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Abstract 

The concentration dependence of the excluded volume effects in polymer solutions is investigated. Through ther- 
modynamic arguments for the intcrpenetration of polymer segments and the free energy change, we show that the 
disappearance of the excluded volume effects should occur at medium concentration. The result is in accord with 
the recent experimental observations. 
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1 Introduction 

Concentration dependence of the excluded volume effects is an unsolved problem in polymer physics. To date, 
only a few theoretical works have been developed on this issue. Quite recently, very reliable experiments with the 
small angle neutron scattering (SANS) have been performed by two research groups [HI!]; they showed that the 
unperturbed coil dimensions arc realized at medium concentration far below the melt state; the unperturbed state 
being retained over wide concentration range. 

In this paper we reexamine the concentration dependence of the excluded volume effects on the basis of the 
classic thermodynamic theory. The starting point of our investigation is motivated by the suggestion of the scaling 
formula by Issacson and Lubensky [3] that a critical concentration may occur for d — 3 somewhere [4] between the 
dilution limit and the non-solvent state above which the excluded volume effects vanish. 

A polymer solution is intrinsically an inhomogeneous system of segment concentration, which is due to the fact 
that monomers are joined by chemical bonds (see Fig. 1). A theory must therefore take this fact into consideration. 
Along this line, first we will make minor amendment of the classic theory of the local free energy in order to apply 
the theory in a more rigorous manner to the disappearance problem. Then, on this basis, we will show that the 
excluded volume effects should really vanish at medium concentration. 

2 Theoretical 

Excluded volume effects [5lll5j of polymer molecules have two different facets: One is the expansion of the polymer 
dimensions and the other is the repulsion between two polymer molecules. These two phenomena, however, can be 
understood by a single thermodynamic property, the osmosis, namely the spontaneous flow of solvent from a more 
dilute region to a more concentrated region. The solvent flowing into from the dilute region will expand polymer 
coils. If all segments are joined by covalent bonds, the osmosis simply leads to the expansion of the polymer coil, 
while if the segments consist of some polymer molecules, the expansion will necessarily lead to the separation 
of those molecules, which will be phenomenologically observed like the core-core repulsion between hard spheres. 
Thus the coil expansion of one molecule and the repulsive interaction between two molecules arc the different 
manifestations of the same phenomenon. 

In order for the osmosis to occur, there must exists concentration fluctuation in the system. Without the 
concentration gradient between the inside and the outside of the coil, no excluded volume effects can occur. Our 



task is then to estimate this gradient and calculate the magnitude of the exclude volume effects as a function of 
polymer concentration. 

In order to investigate whether the transformation of the excluded volume coil to the ideal one occurs at a 
certain concentration, one must begin by evaluating the potential energy change AF of the system under the 
concentration fluctuation. More specifically, one must investigate whether the condition AF -C kT can be fulfilled 
at certain finite concentration. If this inequality is satisfied, the excluded volume effect will be hidden behind 
the thermal fluctuation and never be detected experimentally. Our question is then "To what extent must the 
concentration fluctuation diminish in order for the excluded volume effect to be screened out?" [7l|9]. 

2.1 Local Free Energy: Minor Amendment 

In this section we make minor amendment of the classic theory of the local free energy [5]; the amendment is 
necessary to apply the theory in a more rigorous manner to real polymer solutions. As is well known, the classic 
theory for the local free energy does not take into consideration correctly the interaction between segments on 
the same chain, so that the classic theory is constructed on the basis of the pseudo self-avoiding chains whose 
behavior is virtually equivalent to that of the ideal chain [TTj . This feature has been frequently criticized so far as 
a fundamental deficiency of the Flory theory. 

Let us consider the three dimensional lattice where a site may be occupied by a segment or a solvent molecule. 
Following the standard lattice representation, the random occupation of sites is assumed. Let a site be surrounded 
by z neighboring sites, and let there be 5x segments from different molecules in the volume element SV. The 
number of arrangements of those segments is 

nsv= l[{z^-l){l~f^) (1) 

i=0 

where fi is the probability that a given cell is occupied by the polymer segments when i segments are already put 
in SV^ so that fi = i/Sno where Shq = Sx + Sni, Suq denotes the total number of sites and Sni the number of 
solvent molecules in SV. Zi is a special number introduced in this amendment and is defined by < < z. The 
physical meaning of Zi is as follows: 

There is finite probability that a given segment overlaps with the other segments on the same chain (the 
multioccupation problem). Such unphysical conformations must properly be removed by subtracting from 
the total number Ylii^ ^ 1) of feasible conformations. This is possible, because the number of conformations 
is enumerable in principle. The subtraction can be achieved simply by reducing z to Zi, so that Jli(-^i ~ 1) 
represents the total number of self-avoiding conformations. To date the numerical value of Zi is unfortunately 
unknown |10| . which however is not essential for the present purpose, as is verified below. 

By eq. (1) the local entropy becomes 

5S = SS„iixing + SSQ = k log nsv = k< log (zj + log . °' — — -- > (2) 

[ ^ (Sng^((5no - <5x)! J 

Applying the Stirling formula to the above equation, we have 

{Sx-l ^ 
Sni log vi+Sx - ^ log (zj - 1) ^ (3) 

where SSq = SSqi{Sx = 0) + SSQ2{Sni ~ 0) represents entropy for respective pure components and vi ~ {Suq — 
Sx)/Sno is the volume fraction of solvents. It is clear by eq. (3) that SSoi{Sx = 0) = 0, and SSo2{Sni = 0) = 
—k ^Sx — J2i=o'^ log (zi — Hence we have 

SSmixmg = -k Sui log Vi (4) 

which is exactly the Flory result. It turns out that all the self-avoiding terms arc absorbed into the melting entropy, 
SSo2iSni = 0). Noteworthy is the fact that the mixing entropy SSmixmg of the pseudo self-avoiding chains and 
solvent is exactly equal to that of the genuine self- avoiding chains and solvent; only the standard state must be 
altered as 

pure pseudo self-avoiding chains pure genuine self-avoiding chains. 



This unexpected result does not appear to have been fuhy recognized by theorists up to present. A close investi- 
gation however tells us that the result is by no means surprising, because by eq. (4) SSmixing is a function of the 
solvent fraction alone, but independent of the conformational properties of chains. 

Adding the enthalpy term, SHmixing = kTx6niV2 to the above equation, we have the formula of the local free 
energy in the volume element SV 

5Fmixing = 5H,nixing - T 5 Smiling = kT {log {1 - V2) + XV2} (5) 

where W2 = 1 — f i is the volume fraction of the segments and x the enthalpy parameter. Eq. (5) represents the free 
energy difference between the mixture of the self-avoiding chains and solvent, and the respective pure components. 
Eq. (5) has already been derived by Flory [5]. 

We realize that eq. (5) has deeper generality along with sound physical basis, and hence is applicable equally 
to the excluded volume problem in concentrated systems. 



2.2 Expansion Factor a in Concentrated Solution 

Fundamental force of the coil expansion is the osmosis of solvent 
molecules from a more dilute region, C^aiiey, to a more concen- 
trated region, Chui around the center of gravity of a coil, where C 
denotes segment concentration. In order for the osmosis to occur 
for a given polymer coil, there must exists non-zero free energy 
difference between the inside of the coil [Chui) and the outside 
{Cyaiiey)- If this circumstaucc is realized, positive molecular force 
from the outside region inevitably arises, creating the coil expan- 
sion. 

The size of the expansion factor a is determined by the force 
balance between the osmosis and the retraction force due to rub- 
ber elasticity [5] (see an alternative derivation in Appendix). The 
equilibrium condition is 




Fig. 1: A snapshot of a polymer solution. 



dF/da = dFosmoUc/da + dFeiasUc/da = (6) 
By cq. (5) the free energy of mixing polymer segments and solvent molecules is written in the form: 

AFm ^kT j {log (1 - V2) + x«2} 5r. 



(7) 



where the subscript 1 denotes the solvent, the subscript 2 the solute (polymer segment), and x is the enthalpy 
parameter as mentioned above. Note that 5ni = (1 — V2)5V/Vi with Vi denoting the molecular volume of a solvent 
molecule. Substituting this into eq. (7), we have 



Expanding the logarithmic term, we have 

AFm = — 
Vi 



kT 



(1-W2){l0g(l-1'2)+Xl'2}<5^ 



(1-x) V2 + {l/2-x)vl + -vl + 



5V 



(8) 



(9) 



Since we want to calculate the local free energy difference, AF, between the hill and the valley areas, let us write 
AF in the form: 



AF =AF 
kT 



M,hm 



AF 



AI, valley 



= ^ / i - (1 - X) {-"hill - V^alley) + (1/2 - x) ("m; 



'^valley) + g iyhill ~ ^valley) + 



For W2 ^ 1 which is valid below the medium concentration under discussion, we may write 

kT f 

AF = — y {- (1 - X) {VhUl - V^alley) + (1/2 - x) {vlui ^ ^Llley)} 



SV (10) 



(11) 



Let V2 denote the volume of a polymer segment and we have V2 = V2C. For the Gaussian chain, with the equality 
dV = a'^d{x — a)d{y — b)d{z — c) = dxdydz in mind, we can write generally the segment concentration at a given 
point (x, y, z) in the form 



3 /2 

^-fi{[x-af + {y-bf + [^~cf} (12) 



where N is the number of the segments on a molecule, /3 — 3/2Rg with Rg denoting the radius of gyration of an 
unperturbed chain and (a, b, c) the coordinate of the center of gravity of a polymer molecule; so the summation 
represents the accumulation of segments emanating from a number of different polymers and therefore reflects 
the segment concentration at the point {x,y,z). A graphical representation is given in Fig. 1 which shows that 
a polymer solution is the typical inhomogeneous system. In contrast to the homogeneous system of a monomer 
solution, there is wild fluctuation AC = Chui — Cyaiiey in the solution. AC is a function of a; it rapidly decreases 
with increasing a, because segments pervade more deeply the whole system as a increases. 
It is useful to extract the pre-factor from eq. (12) and recast it in the form 

13 



Cix,y,z)^N[^\ G{x,y,z) (12') 
so that the quantity 

G(x,y,z)= e-''«--'')'+(^'^)'+(^-=)'> (13) 

{a,6,c} 

is a function independent of a. Now cq. (11) may be rewritten in the form 

y2 fffi , ^fPV^^rr. ^ N AW /o ^ /3 ^ ^ 



AF = NkT^^ jjj)^-(i-^)^t_j ^Gmu - Gvalley) + V2N (1/2 -X)y-^j {Glu - Gl^aey) I dxdydz 

(14) 

The first term of cq. (14) is, by eq. (13), independent of a, so it vanishes by the differentiation with respect to a. 
We have then 

V2' _Jfi'^' 



d^FosraoUc/da = -3N^kT-2^^ (1/2 - X) J j JJJ {Glu - G^iey) dxdydz (15) 

Note that for any fluctuation model of polymer solutions, we can establish one-to-one correspondence between the 
Ghill and the Gvaiiey by properly setting the defined spaces. Thus it is sufficient to take account of the elastic 
force per one molecule (this point will be discussed in Sec. 4). The classic work [T^ showed that dAFf.iastic/da = 
3kT [a — 1/a). The force balance between the osmosis and the rubber elasticity is obtained by substituting these 
results in eq. (6), namely 

37V2fcT^ (1/2 - X) (0 /// {GLi - GLuey) dxdydz = 3kT{a - 1/a) (16) 
By rearrangement, we arrive at the expression (see an alternative derivation in Appendix): 

(1/2 - X) (0 /// {GLi - GLuey) dxdydz (17) 

For the dilution limit, we have Ghui — > e'^^^^'^^^'^''^^ and G^aiiey — > and eq. (17) exactly reduces to the Flory 
result. 

There are important implications in eq. (17): (i) A polymer coil is ideal (a = 1) at the theta point (x = 1/2) 
as is well established already, (ii) for a large Vi limit, the coil is nearly ideal, and (iii) the fifth power rule of a is 
still valid, but (iv) the coil must also be ideal at the point of Ghui = Gyaiiey in which the concentration fluctuation 
disappears. 

2.3 Excluded Volume u in Concentrated Solution 



Prior to the derivation of the excluded volume u in the concentrated solution, let us review briefly the classic work. 
Let a system contain only two polymer molecules. Consider the interpenetration of the two polymer molecules. 



Let a distance be L between their centers of gravity. The chemical potential change (Ai^ = Ai7 — TAS") bringing 
these molecules close from L = oo to i is 

^F{L) - 2kT (1/2 - x) {y^lVi) J PkPe dV (18) 

where pk and pi are the respective segment density of the two polymer molecules k and £ within the small volume 
element SV, which has the form: 

^ N ( \-fi-'-+yi+^i) (19) 



Then the excluded volume u can be calculated by the equation 

'47rL2fi_e-Ai=^(L)AT\^^ (20) 



The parameter x is, according to the definition, positive for poor solvents but negative for good solvents. It follows 
that the interpenetration of the two polymer coils is strongly hindered in good solvents resulting in the molecular 
repulsion. We derive in the following the corresponding free energy change in the concentrated polymer solution. 

2.3.1 Free Energy Change 

Suppose that a single polymer molecule k is added to the concentrated polymer solution. There are two cases for 
the molecule to be located in the solution: one is the case where the molecule is put within the Chui region and the 
other is the case where it is put within the Cvaiiey region. The free energy difference between these two locations 
will correspond to the free energy difference between overlapping state and non- overlapping state in concentrated 
solutions. Now we can generalize the Flory excluded volume theory to include the concentrated system. The 
quantity in the small volume element 5V is 

5{AFk.c) = S{AFu^Cum)~ S{AFk,c.^„J = kT{SV/Vi) { (1 - pkV2 - ChuiV2)log{l - pkV2 - C^/Fs) 

- (1 - ChUlV2) log (1 - ChUlV2) - (1 - PkV2 - G,aUeyV2) log (1 - pkV2 - C.,aUeyV2) 

+ (1 - C.aUeyV2) log (1 - C^aUeyV2) ~ ^XPkVi {CmU ~ Cvalley) } (21) 

Eq. (21) can be written in the series form: 

5{AFk,c) = 2A:T{(l/2 - x) + O} pu {Cmu - C^aiiey) {V^ /Vi) 5V (22) 

where the symbol O represents the higher terms of the series and a function of Chui, Cyaiiey and pk- Below the 
medium concentration, these terms are negligible and eq. (22) reduces to 

(5(AFfc,c) = 2kT (1/2 - x) Pk (Chui - C^aiiey) {Vi/Vi) 6V (23) 

Let L be a distance between the centers of pk and Chui ■ Then Cyaiiey is also a function of L, because the distance 
between Chui and Cyaiiey is fixed, on average, for the solution of a given concentration. Applying the equality 
SV = a^dxdydz to eq. (23) and integrating over the defined spaces, we have the free energy difference as a function 
of L 

AFk^ciL) = 2kT (1/2 - x) [V^ /Vi) a"" j j j pu [Chai - Cyaiiey) dxdydz (24) 

In the limit of the infinite dilution, we have Chui pe and Cyaiiey ^ 0, and we recover the classic equation (18). 
Using the Gaussian approximation, we may recast the above equation in the form: 

AFk,c{L) = 2N^kT (1/2 ~ x) {V^ /Vi) 1^-^ jjj e-^(^"+y"+^i) {Gmu - G valley) dxdydz (25) 



ill 



where G = J2{ab c}^~ +(^f-c) }^ ^^s the concentration fluctuation decreases, namely, Gh 

Gvaiiey 0, then AFk^c ^ 0. And it follows from cq. (20) that u 0. The excluded volume disappears in 
parallel with the disappearance of the fluctuation. 



3 Solution According to Lattice Model 



3.1 Evaluation of Fluctuation 

By the equations (17) and (25), we have a general statement that the excluded volume effects are a strong lunction 
of the concentration fluctuation. This is a very natural conclusion, because the osmosis can occur only in the 
presence of the concentration gradient. Our task is then to evaluate the local concentration gradient around a 
given polymer coil. 

Eqs. (17) and (25) are formal solutions not easy to solve, since the relative coordinates (a, b, c) of individual 
molecules can not be specified in real solutions; moreover the boundary conditions of the integral terms in eqs. 
(17) and (25) are not clear. In this sense, no rigorous calculation seems possible. However it is possible to extract 
essential features of the equations by making use of the lattice model. In this paper we show the solution for eq. 
(17) only, since eq. (25) is more complicated and requires heavy calculation that seems to exceed the ability of 
computers available. 




Fig. 2: Representation of the simple cubic lattice: Edges Fig. 3: Representation of the concentration fluctuation 

are occupied by polymer molecules having the Gaussian in the solution of Gaussian polymers arranged on the 

distribution. The C maximum and the minimum points simple cubic lattice. A peak represents the Chui and a 

lie on the z = x plane. bottom the Cvaiiey'- the center of a peak corresponds to 

the center of gravity of a chain. 



Consider the simple cubic lattice with the unit lengths {p x p x p). Polymer molecules are arranged on every 
sites with the Gaussian segment distribution. The relative coordinates of the molecules have then integer values 
of the form: (a, 6, c) = {ip, jp, kp) {i,j,k = — oo, • • • ,— 1,0,1,-- - ,+oo). The segment number density (/A ) at a 
given coordinate {x,y, z) can be calculated by the simple summation of equal steps: 



+ 00 . o X 3/2 



+ 00 ^ ^ ^ 3/2 

{i^j.k) — ~oo 
^ X 3/2 



=Nl^\ G{x,y,z) (26) 
To find the maximum and the minimum points, differentiate eq. (26) to yield 

Since y and z are independent of each other, the solution for the equality dC = must satisfy 



dC\ fdC\ fdC 



Here 



^ ) = const ^ (x - ip) e-/5{(-'rt^+(y-^rt^+(-M^} ^ q (29) 



dx 



+ 00 



= (28) 



Obviously x = ip satisfies eq. (29). The other solutions are x = {l/2 + £)p (£ = — cx), • ■ • , —1,0, 1, ■ • ■ , +oo). From 
these, we have 



C,nax = {x = ip, y = jp, z = kp) 

C™„ = {x = (1/2 + i)p,y = (1/2 + m) p, z = (1/2 + n) p) (30) 

where £, m, n = —oo, • • • , —1, 0, 1, • • • , +oo. Thus Cmax and Cmin lie on z = x + np planes. In Fig. 2, an example 
{z = X plane) of those planes is illustrated. The mean number density of segments is directly calculated by the 
equation: 

C = ^ i/k') (31) 

so that the mean polymer volume fraction cf) is 

- N , , 

<i>=—V2 32 

An important quantity to be evaluated is the integral term appearing in the final expressions (17); this measures 
the 3-dimensional density fluctuation in the system. We put 

^ = /// {Glm ~ Glaiiey) dxdydz (33) 

As can be seen from Figs. 1-3, the Ghui and the Gyaiiey areas are discontinuous; its domains can be partitioned 
into a number of intervals. By the lattice symmetry, we may define the intervals as [— p/4,p/4] for Ghui for each 
axis and [p/4, 3p/4] for G valley This choice is simply a matter of convenience for numerical calculation. Then the 
above integral may be specified as 

r r rP/4: r r r3p/i 

Ja^ G^ dxdydz - / / / G^ dxdydz (34) 

JJJ-p/4 JJJp/i 

We show in Fig. 4-a the concentration fluctuation on the x = y = z line that goes through the centers of the 
Chill and Cyaiiey The curves were calculated according to eq. (26) for a = 1 as a function of the mean polymer 
volume fraction 0, modeling polymethyl methacrylate (PMMA: N = 1000) solutions. We have assumed the size 
of the segment to be equal to that of the repeating unit so that V2 = 140 A. The light-brown peak with = 
represents the isolated polymer molecule in the dilution limit. As one can see, the concentration fluctuation decays 
rapidly with increasing (j) and vanishes for ^ > 0.2. 

In Fig. 4-b, the quantity Ja/N is plotted as a function of 4> and N. The dotted line shows the result for 
N = 1000 and the solid line the result for N = 5000. The function Ja/N also decreases strongly with increasing 
concentration and vanishes at medium concentration. 
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Fig. 4: Concentration fluctuation in polymethyl methacrylate solutions (V2 
X = y — z line (see also Fig. 2) as a function of the volume fraction <j) of the polymer (A'" — 1000). (b) Numerical 
solution of Ja/N as a function of (j>; Ja was calculated according to eq. (34) together with eq. (32); dotted line (• • •): 
N = 1000 and solid line (-): N = 5000. 
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3.2 Expansion Factor a as a Function of Concentration 



With the help of the estimation of Jq, we can now solve eq. (17) as a function of the segment concentration. 
Employed parameters are listed in Table 1. The enthalpy parameter x is hard to estimate. Here we use the values 
of X = 0.3 for PMMA— chloroform and x = 0-2 for PE— 7i-nonadecane which is an assumed value at 150 °C by the 
extrapolation of the HandBook data [17] . 

The calculation results are illustrated in Figs. 5 and 6 to be compared with the observed points by Cheng, 
Graessley and Melnichenko jT| and those by Westermann, Willner, Richter and Fetters [2]. As one can see, the 
expansion factor a drops strongly with increasing polymer concentration and falls to the unperturbed value at the 
medium concentration in accordance with the experimental observations [HE] ; the results are also consistent with 
the prediction of the Issacson-Lubensky scaling formula [Sill]. 

According to the present calculation (Figs. 5 and 6), it is found that the essential features of the two systems 
of PMMA and PE are very alike: both of the systems show the swollen-to-unperturbed coil transition at the 
medium concentration. Comparing the two systems, it is found also that the coil expansion is less pronounced 
in the PE— n-nonadecane system than in the PMMA— chloroform system. The reason can be found by inspecting 
eq. (17); it simply comes from the special combination of the polymer and the solvent, namely n-nonadecane with 
the large molecular volume Vi and polyethylene with the small segment volume V2- According to eq. (17), such a 
combination necessarily depresses the expansion factor a to a lower level. 



Table 1: Basic parameters of polymer solutions [l][2] 





parametersj_^ 


notations 


values 


polymethyl methacrylate (PMMA) 


volume of a solvent (CHCI3) 


Vi 


134 A' 




volume of a segment (C5O2H8) 


V2 


140 A' 




degree of polymerization 


N 


5900 




Flory characteristic ratio 


C p 


9.2 




mean bond length 


I 


1.56A 




enthalpy parameter (25 °C) 


X 


0.3 


polyethylene (PE) 


volume of a solvent (n-Ci9H4o) 




569 A' 




volume of a segment (C2H4) 


V2 


49 A' 




degree of polymerization 


N 


1000 




Flory characteristic ratio 




7.7 




mean bond length 


I 


I.54A 




enthalpy parameter (150 °C) 


X 


0.2 



4 Discussion 

The lattice model is a strong approximation where the dynamic aspect of real solutions is completely neglected. 
However, there are two cases in which this static model is expected to be a good representation of real solutions: 
(i) the lattice model will predict correctly the behavior in dilute solutions where polymer coils are, on average, far 
apart from each other so that the interaction among chains is weak; (ii) the model will predict correctly the location 
of the swollen-to-unpcrturbed coil transition at which the concentration fluctuation disappears. The agreement 
with the recent observations supports this reasoning. 

One of the important questions throughout the present work was whether one can put the equality between the 
l.h.s and the r.h.s. terms in eq. (16); i.e., whether one can equate the osmotic pressure due to the density gradient 

t'^ The unperturbed size is calculated by the equation: = ^CpN£^{?, where ^ is the bond number per one repeating unit. The 
numerical estimation of Rg is in accord with the observed values 1181 . 
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(p (volume fraction of polymer) 

Fig. 5: Expansion factor vs <t> plot for PMMA-CHCI3. 
Solid line (— ): theoretical line by eq. (17) for x = 0-3- 
Dotted line (•■•)• theoretical line by eq. (17) for x — 0- 
Open circles (o): observed points by Cheng, Graessley 
and Melnichenko. 
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Fig. 6: Expansion factor vs cj> plot for PE— n-Ci9H4o. 
Solid line (— ): theoretical line by eq. (17) for x = 0.2. 
Open circles (o): observed points by Westermann, Will- 
ner, Richter, and Fetters. 



exactly with the retractive force of one whole chain. To examine this problem, we must take the following fact into 
consideration: the osmotic phenomena within polymer coils are very different from the ordinary osmosis |19| . The 
ordinary osmosis occurs between two systems separated by the semipermeable membrane, in which solvent molecules 
diffuse from a dilute system to a concentrated system, thus lowering the solute concentration in the concentrated 
system, whereas raising that in the dilute system. This process of the change of the solute concentrations on each 
side of the membrane takes place as a result of the volume change of the respective systems. 

In contrast, there exists no semipermeable membrane in polymer solutions so that no volume change of the 
respective systems can occur. The change of the concentrations in the two regions (dilute and concentrated) can 
occur only through the interchange of the space coordinates between polymer segments and solvent molecules, 
which leads to the coil expansion. The expansion of the polymer dimensions is therefore similar to the dissolution 
of a solid material into a solvent. The coil expansion resulting from the interchange of the coordinates tends to 
reduce the density fluctuation in the system, whereas the coil contraction will augment it. Only one way for the 
polymer solution to reduce its own fluctuation is to expand the polymer dimensions. Such adaptation of the coil 
dimensions should occur uniformly over the whole system because of the absence of the semipermeable membrane, 
thus validating the equality (16). Note that in the alternative derivation shown in Appendix, this problem on 
the equality (16) does not arise explicitly because of the delta function approximation for the segment density 
distribution. 

There exists another question, "what do we mean by the unperturbed (or ideal) chain?" To answer this question, 
let us return to the basic assumption of the present work. We have calculated the free energy difference AF by 
subtracting the free energy in the dilute region from that in the concentrated region. Thus the term "unperturbed" 
necessarily refers to the configuration under the condition that satisfies Cmax = Cmin- By the result of eq. (30) 
it corresponds to the configuration at p = 0, namely the configuration at infinite concentration (C — > 00) where 
all excluded volume effects should vanish rigorously as discussed earlier [3] . Hence the term "unperturbed" defines 
a standard state for the configuration of a chain composed of mathematical dots and lines (with no volume and 
no thickness) and immersed in the true athermal solvent. In that hypothetical limit, the segment should take the 
distribution of the form [20] : 




where cuk = t^. + (1 — tk)^, t = k/N and R is the end-to-end distance and ~ 6Rg^. Usually this distribution 
is approximated by the corresponding Gaussian distribution as we have done in the text. Note that the above 
definition for the unperturbed chain is clearly different from that used in the classic theory [5] where the standard 
state has been taken to be a pure polymer (melt state of self-avoiding chains). The two definitions are, however, 
virtually the same, since a chain obeys in effect the random flight statistics over all concentration range above the 
medium concentration, as we have seen in Figs. 5 and 6. 
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Conclusion 



1. We have made a minor amendment for the local free energy, showing that the classic theory has deep generality 
and sound physical basis. 

2. Making use of the thermodynamic arguments, wc have developed the extended theory of the excluded volume 
effects, showing that the expansion factor a is a strong function of the concentration fluctuation eq. (17). 

3. According to the solution of the lattice model, the fluctuation is maximum in the dilution limit, but decreases 
strongly as the polymer concentration increases, and vanishes at the medium concentration Fig. 4. 

4. In parallel with the behavior of the fluctuation, the excluded volume effects manifest themselves most pro- 
nouncedly in the dilution limit, but decay rapidly with increasing concentration and vanish at the medium 
concentration. 

The theoretical results are in good accord with the recent experimental observations [Tl[2] Figs. 5 and 6. 



6 Appendix 

An Alternative Derivation of Eq. (17) 



Let W{s) be the unperturbed segment distribution around the center of gravity. We introduce the partition function 
for the excluded volume chain l6l[2T]: 



W{s) cxp (^--^^ Airs^ds 



(Al) 



where V(s) represents a potential function. Then the perturbed segment distribution can be formulated as 



P{s) = Iwis) exp (-^ ] Ans^ (A2) 



Z ' ' V kT 
Making use of the Gaussian approximation for W{s), we have 

Rijp{s)ds je^^^j-t^-my-dt 

In eq. (A3), we have made the variable transformation: ~ s^/R^. Independently we have another equality 



= ^ = = -r ^^rv< (A3) 



a 



t6{t~a)dt 

a' = (A4) 

/ -6{t-a)dt 

where d{t — a) signifies the delta function peaked at i = a. If we identify eq. (A3) with eq. (A4), this amounts 
to making the approximation: S{t — a) « exp ^— f^^ — t^- Our remaining task is then only to evaluate the 

maximum point of the function, exp ( — f^^ — ^Vr') so that 



^^log|exp(--t 



(A5) 



Here we use AF{a — > t) in eq. (14) as the potential function V{t). After some rearragement, we obtain 

/(l/2-x)f^'' 

Vl V TT 



a'-a^= N^Yl_ (1/2 - x) ( § ) / / / (GL; " G^aey) dxdydz (A6) 



which is just eq. (17) in the text. 
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